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Abstract 

The physical properties of the commonly used second-order closure models 
are examined theoretically for rotating turbulent flows. Comparisons are made 
with results which are a rigorous consequence of the Navier-Stokes equations 
for the problem of fully-developed turbulent channel flow in a rapidly 
rotating framework. It is demonstrated that all existing second-order 
closures yield spurious physical results for this test problem of rotating 
channel flow. In fact, the results obtained are shown to be substantially 
more unphysical than those obtained from the simpler K-e and K-2, models. 
Modifications in the basic structure of these second-order closure models are 
proposed which can alleviate this problem. 


Research was supported by the National Aeronautics and Space 
Administration under NASA Contract No. NAS1-17070 while the author was in 
residence at the Institute for Computer Applications in Science and 
Engineering, NASA Langley Research Center, Hampton, VA 23665-5225. 
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1. Introduction 

The ability to accurately predict rotating turbulent flows could have 

a wide variety of important applications ranging from the analysis of 

turbomachinery to the description of turbulence in the atmosphere or in 

oceans. However, with the exception of some recent work on the limiting 
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case of two-dimensional turbulence , there appears to have been little 

work which accounts for the effects of rigid body rotations on turbulence 

modeling by a direct analysis of the Navier-Stokes equations. In fact, 

much of the work in modeling rotating flows has been conducted utilizing 

the K-e (or K-il) model of turbulence which requires the effects of rotations 

to be accounted for by the use of a variety of ad hoc empiricisms (c.f. , 

Majumdar, Pratap, and Spalding , Howard, Patankar, and Bordynuik , and 

Galmes and Lakshminarayana^) . While such approaches can be useful in 

correlating experimental results for a particular class of rotating turbulent 

flows, they do not form the basis for a general theory which is needed if 

models are to be developed that truly have a predictive value. 

The purpose of the present paper is to examine the consistency of 

various second-order closure models with results which are a rigorous 

consequence of the Navier-Stokes equations in a rapidly rotating framework. 

Although there have been a few previous studies of rotating turbulent flows 
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using second-order closure models (c.f., So , So and Peskin , and Mellor 

9 

and Yamada ), only a narrow range of flows have been considered so that 
no definitive conclusions could be drawn about the correctness of the models. 
In this paper the test problem of fully-developed turbulent channel flow 
in a rapidly rotating framework will be considered. It will be proven, 
as a rigorous consequence of the Navier-Stokes equations, that the Reynolds 
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stress tensor for this problem must be two-dimensional (as a direct conse- 
quence of the Taylor-Proudman theorem"^) and must have a non-zero Reynolds 
shear stress in the plane of the flow. The commonly used second-order 
closure models yield completely opposite results (i.e., a three-dimensional 
Reynolds stress tensor with vanishing Reynolds shear stresses), and, hence, 
are fundamentally inconsistent with the Navier-Stokes equations for turbulent 
channel flow in a rapidly rotating framework. It is also demonstrated 
that the results obtained from these second-order closures are substantially 
more unphysical than those obtained from the simpler K-e or K-£ models. 
Modifications in existing second-order closure models which can alleviate 
this problem will be explored. By the addition of one term to the second- 
order closure obtained by Haworth and Pope^ from a generalized Langevin 
model, it will be shown that consistency with the Navier-Stokes equations 
in a rapidly rotating frame can be achieved. Other possible modifications 
will be discussed briefly in the last section along with the prospects for 
future research. 


t 
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2. Second-Order Closure Models and Rotating Channel Flow 

The turbulent flow of a homogeneous and incompressible viscous fluid 
in a rotating frame of reference will be considered. The velocity field v 
and pressure P will be decomposed into ensemble mean and fluctuating 
parts, respectively, as follows 

v = v + u, P = F + p (1) 

Of course, the decompositions in Eq. (1) are solutions of the mean and 

fluctuating parts of the Navier-Stokes equations and continuity equation 

12 

which, in a steadily rotating frame of reference, take the form 

3v _ _ _ 2— — 

+ vVv = - VP + vV v - V*T - 2Q x v (2) 

at ~ - ~ ^ ~ ~ 


3 u _ _ 2 

■7-— + v*Vu = - u*Vu - u*Vv - Vp + vV u + V* x - 2Q x u (3) 

dt ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ^ 


V*v = 0 (4) 

V-u = 0 (5) 


where v is the kinematic viscosity of the fluid, is the angular velocity 
of the reference frame relative to an inertial framing, and t is the kinematic 
Reynolds stress tensor given by 


x = uu (6) 

It should be noted that "p in Eq. (2) is the modified mean pressure which 
includes the centrifugal and gravitational body force potentials. 
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The Reynolds stress transport equation is obtained by taking the 

ensemble mean of the symmetric part of the outer product of Eq. (3) with u. 

12 

This equation takes the form 


DT 


k£ 


3v, 


Dt 


+ T. 


km 9x 


+ T 


m 


£m 3x 


m 


3x 


5- 'W(> 


m 


- (u 


9 p 
k 9x, 


+ u 


9u i 9u o o 

l ^ - 2v ST ST * vV T U 


m 


m 


+ 2 


£ , T „ + 
mkn m nJc 


2e .fix. 
mien m nk 


(7) 


in a rotating frame of reference. In (7), Cartesian tensor notation is utilized 
where the Einstein summation convention applies to repeated indices and 
represents the permutation tensor. Second-order closure models are obtained 
when closure relations are provided that tie the higher-order turbulence 
correlations in (7) to the Reynolds stresses (along with their spatial 
gradients), the mean velocity gradients, and the length scale of turbulence. 

The closure relations used for the higher-order correlations (i.e., the 
triple velocity correlation, pressure gradient-velocity correlation, and the 
dissipation rate correlation) in an inertial framing are of the general form 


Wm = 


( 8 ) 


U k 9x + U £ 9x " C l (T km 9x + T £m 9x 5 + 
36 k mm 


(9) 


9u k 9u £ 

2V 9T fcT = 
m m 


( 10 ) 
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13 14 

respectively (see Mellor and Herring , Launder, Reece, and Rodi , and 
Lumley^) . In (7)- (10), is a dimensionless constant, SL is the length 
scale of turbulence which must be obtained from a separate transport equation, 
and 


D 

Dt 


3_ 

at 


+ v*V 


D 


i. 3v k 

- — + 


k£ 2^3x 


£ 


3x, 


( 11 ) 

( 12 ) 


are, respectively, the mean substantial derivative and rate of strain tensor. 

2 

By extending (8)- (10) to rotating frames of reference (see Speziale ) and then 
substituting the results into (7), it is a simple matter to show that all 
existing second-order closure models are of the form 


DT^ + ^ c i> sr * ni 

m m 


)“¥■ W 1-lZ-V 

m 


- A kil (T,VT,p,£) - B k£ (T,D,£) + W 2 T kii 

+ (2-C_ ) (e . Q t n + c n ft t . ) 
v 1 mkn m ni m£n m nk' 


(13) 


in a rotating frame of reference. It should be noted that in the Rotta- 
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Kolmogorov model the constant = 0 and in the Launder, Reece, and Rodi 

and Lumley^ models, is a non-zero constant which does not equal 2. 

Hence, the Coriolis term on the right-hand-side of (13) survives in all 


existing second-order closures where it constitutes the only non-inertial term. 
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Now, we will consider the test problem of fully-developed turbulent 
channel flow in a rapidly rotating framework (see Figure 1). This fully- 
developed channel flow is maintained by a constant axial mean pressure gradient 


9P 

8x 


- G 


(14) 


while the channel is subjected to a steady spanwise rotation with the 
angular velocity (see Figure 1) 


n = ftk 

The mean velocity field is given by 
v = U(y)i 

and the Reynolds stress tensor is of the form 


(15) 


(16) 


T = T(y) 


(17) 


since the flow conditions are fully-developed. Hence, the axial component 
of the Reynolds equation (2) in a rotating channel flow is given by 


v 




+ G 


0 


(18) 


Since it is well known that the axial velocity profile in fully-developed 
turbulent channel flow is not parabolic like its laminar counterpart, it is 
clear that 


T 

xy 


t 0 


(19) 
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(also see the experimental results of Johnston, Halleen, and Lezius^). 
The y and z- components of the Reynolds equation (2) take the form 


-XT dx 
9P _ yy 

9y dy 


- 2 ftU = 0 


( 20 ) 


3z dy u 


( 21 ) 


Since there is no forcing in the z-direction. 


I=° 


and (21) can be integrated to yield the result 


x =0 
yz 


( 22 ) 


since x must satisfy the no slip condition at the channel walls. Similarly, 
it can be shown that 


x =0 
xz 


(23) 


The vorticity transport equation in a steadily rotating framework (where 
co = V xv) takes the form ^ 


Deo 2 

~ = co*Vv + vV co + 2ft *Vv 
Dt ~ ~ ~ ~ ~ ~ — - 


(24) 


A rapidly rotating framework is defined as one in which 


ft t » 1 

o 


( 25 ) 
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where t Q is the time scale of the turbulent fluctuations. For a rotating 
channel flow where = ftk, (24) can be written in the equivalent form 


Deo 

Dt 


to* Vv 


v „2 

a v “ 


+ 2k*Vv 


(26) 


As a result of (25), equation (24) yields the approximate constraint (in 
a rapidly rotating framework) that 


2k* Vv = 0 (27) 

It, thus, follows that 
3v 

"9z = ° (28) 

Since the walls of the channel are parallel to the axis of rotation, (28) 
has the exact solution 

v = v(x,y,t) (29) 

which is valid throughout the channel . Hence, the flow must become two- 
dimensional (in a strong approximate sense) as would be expected from the 
Taylor-Proudman theorem. 10 Since, v = U(y)i, it is clear that the fluctuating 
velocity is of the approximate form 

u = u(x,y,t) ( 30 ) 

in a rapidly rotating channel flow. It will now be shown that since there 
is no mean forcing along the axis of rotation (i.e., since = 8P/3z = 0) it 
follows that 


u 


0 


(31) 
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This result can be easily seen by examining the z-component of the fluctuating 
momentum equation (3) which reduces to the form 


3u 3u 

z 77 z 

+ U -r 

3x 


at 


3u 3u 0 

z z 2 

U tt U -TT — + vV u 

x 3x y 3y z 


(32) 


since 3p/3z = 0 (the flow is unforced in the z-direction) . Once U, u^, and 
u are determined from the x and y-components of the momentum equations (2) 
and (3) , equation (32) constitutes a linear partial differential equation 
for the determination of u^. This equation has the simple solution u^ = 0 
when solved subject to the no slip condition at the channel walls. Hence, 
the solution 


T 

XZ 


u u 
X z 


= 0 


(33) 


T 

ZZ 


u u 
z z 


0 


(34) 


is consistent with the Navier-Stokes equations. It is thus clear that, 

as a rigorous consequence of the Navier-Stokes equations, the Reynolds 

stress tensor in a rapidly rotating channel flow takes the two-dimensional form 



t (y) 

XX J 

t (y) 

xy 

0 


T = 

x (y) 
xy 

X (y) 

yy 

0 

(35) 


_ 0 

0 

0 



This differs from T for turbulent channel flow in an inertial framework only 

in that t is zero -- a result which is a direct consequence of the Taylor- 
zz 

Proudman reorganization that occurs in a rapidly rotating framework. 
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It will now be demonstrated that all existing second-order closure 
models yield results that are fundamentally inconsistent with (35) and, hence 
are in serious violation of the Navier-Stokes equations for rotating channel 
flow. By establishing a one-to-one correspondence between (x,y,z) and 
(XpX 2 >X 2 ), the general form of second-order closures given by (13) reduces 
to the equation 


1 [ D V . . 9v <! 3v 

a [ Dt- + '‘-'ll'VF* T Hm to 

m m 


k,, 1 

JJ 3x 

m 


- 5 A kik - h v 2 T. 


ft 


k£ 




(36) 


where we have made use of the fact that ft = ftk and we have divided through 
by ft. Taking the limit of (36) as ft -»- we conclude that all existing 
second-order closure models yield the constraint 


( 2_C l) ( £ 3kn T n£ + SjlnW " 0 


(37) 


for a rapidly rotating channel flow. As noted earlier, 

2 - Cj / 0 


(38) 


for all of the existing second-order closures (in fact, Cj = 0 for the 
Rotta-Komogorov model and Cj < 1 in the Launder models 14 ). It is then a 
simple matter to show that (37) has the general solution 


T = T 
xx yy 


xy 


(39) 


0 


(40) 
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Hence, all existing second-order closures predict a state of transverse isotropy 
for turbulent channel flow in a rapidly rotating framework --a result which 
is in conflict with the Navier-Stokes equations as demonstrated earlier. 

In fact, this constitutes a completely spurious physical result since it 
is well known that a turbulent shear flow must be accompanied by a non- zero 
turbulent shear stress. Furthermore, since the Coriolis term vanishes in 
the zz-component of all existing second-order closures (see Eq. (13)), it 
is clear that these models predict that 

T zz i 0 (41) 

as they do for a turbulent channel flow in an inertial framing. This result 

is inconsistent with the Taylor-Proudman theorem for rotating channel flow 

as discussed earlier. These inconsistencies arise because all existing 

second-order closure models violate the principle of material frame-indifference 

in the limit of two-dimensional turbulence -- a result which is a rigorous 

1 17 

consequence of the Navier-Stokes equations as proven by Speziale. ’ To 
be more specific, the inertial terms in (13) do not vanish in a two- 
dimensional turbulence unless it is isotropic. In the next section, it will 
be shown how these deficiencies can be remedied. 
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3. Alternative Turbulence Models For Rotating Flows 

As a prelude to examining alterations in existing second-order closure 
models which can remedy the deficiencies discussed in the previous section, 
we will consider the consistency of the simpler K-e or K-5. models of 
turbulence with the Navier-Stokes equations for rotating channel flow. In 
the K-e or K-£ models of turbulence, the Reynolds stress tensor takes the 
form 


T k£ = 3 K 6 k£ 


i 9v k 
aK 2 £(-r— + 
^ 9x„ 


3v^ 

9x, 


(42) 


where 


K = *“ T 
2 mm 


(43) 


is the turbulent kinetic energy per unit mass which is obtained from a 
modeled version of its transport equation, % is the length scale of 
turbulence, and a is a dimensionless constant. For the K-£ model of turbulence, 
l is either specified algebraically or is determined from a modeled 
transport equation obtained from the contracted form of the evolution 
equation for the two-point velocity correlations. In the K-e model, 
the length scale is given by 

* 3/2 

1 = ~r~ (44) 

where e is the dissipation rate of the turbulence which is obtained from a 
modeled version of its transport equation. 14 Since the tensorial dependence 
of ( 42) is only on the Kronecker delta and the symmetric part of the mean 
velocity gradients which are frame- indifferent, these traditional K-e 
or K-£ models are of the same form in a rotating frame as in an inertial 
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framing. The only way that the effects of rotation are accounted for is 

through a change in the length scale which is usually done with ad hoc 

empiricisms (see Howard, Patankar, and Bordynuik 5 ) . This result (along 

with the fact that the mean flow is unidirectional so that the Coriolis 

17 

acceleration in (2) affects only the mean pressure ) guarantees that 
the basic structure of the Reynolds stress tensor predicted by the K-e or 
K-& models will be the same in an inertial or a rotating frame of reference. 
To be more specific, the Reynolds stress tensor obtained from K-e or K-£ 
models of turbulence will be of the form 


X T 

xx xy 

X T 

xy yy 

0 0 


0 

0 


T 

ZZJ 


( 45 ) 


(where x = x = x ) for rotating channel flow or for channel flow in 

XX jJ Z Z 

an inertial framing. Of course, the specific values of the non-vanishing 

components of x can be different in the rotating frame as a result of changes 

in the length scale. While (45) is not correct for a turbulent channel flow 

in a rapidly rotating framework (since x^ f 0) , it is substantially less 

unphysical than results obtained from existing second-order closures where 

x is zero, 
xy 

More recently, Galmes and Lakshminarayana^ developed an algebraic 
model for rotating flows which is given by 


2 2 K 1 

x. 0 = -=• K<5, „ - 4—— (1-4 C 0 )(e, n x „ + e n fix,) 
k£ 3 k£ Cj£ 2 2 A kmn m n£ £mn m nk J 


Kcwy , 


3v 


T T „ 


m 


3v, 


3v, 


- (T. 


(he L 3 l mn 8x u k& *■ L km Bx + T £m 9x 
1 n mm 


)] 


( 46 ) 
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where and are dimensionless constants that take on the approximate 
values of 1.5 and 0.6, respectively. For rotating channel flow where 

fl = Ok 

Eq. (46) can be written in the equivalent form 


!2 x k)T 3 a ! kl * C £ ^ ' 2 C 2 Jf e 3kn x n£ + c 3£n x nk J 

K0^ , jV, 

fiC, e 3 mn 3x k£ ^ km 3x £m 3x ^ 
1 n mm 


which, in a rapidly rotating framework (i.e., for £}-»-«»), yields the approximate 
constraint 


£ 3kn X n£ + e 3£n X nk 


(48) 


Eq. (48) requires that t ■ = x and t =0. Hence, the algebraic Reynolds 

xx yy xy 

stress model for rotating flows proposed by Galmes and Lakshminarayana^ 
yields the same spurious result for a rapidly rotating channel flow (i.e., 
a vanishing Reynolds shear stress in a turbulent shear flow) as that predicted 
by existing second-order closure models. This problem arises since the 
inertial term in (47) does not vanish for a two-dimensional turbulence 
and, hence, (47) is in violation of the principle of material frame-indifference 
which is a rigorous consequence of the Navier-Stokes equations in the two- 


dimensional limit. 
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It will now be demonstrated that these problems can be overcome by making 
a small modification in the second-order closure recently obtained by 
Haworth and Pope 11 from a generalized Langevin model. In this Langevin 
model, a Reynolds stress closure is obtained which, in an inertial framing, 
takes the form 


Dx 

Dt 


k£ 


+ 


T km 3x + T £m 3x 
m m 


■x — ( u u.uj 
3x m k £ 


m 


+ G . t « + G 0 x , + C e6. „ 
km m£ £m mk o k£ 

where C q is a dimensionless constant and G is a second rank tensor which 
is assumed to be of the general functional form 


(49) 


G — G(t, Vv, e) 


(50) 


It should be noted that in this Langevin model it is not necessary to model 
the triple velocity correlation uuu; this term is closed since it arises 
from the convection in physical space of the joint probability density 
function which is obtained from a separate transport equation. Haworth 
and Pope 11 arrived at a model for G of the form 

! 3v 

G k£ “ (2+e i ) 377 " 2K B 1 dT T m£ +H k£ ( 7’5’ e;) (51) 

£ m 

where 3^ is a dimensionless constant and H is a frame- indifferent function 
since it depends only on the frame-indifferent quantities x, D, and e. Hence 
the Langevin model (49) is, in an inertial framing, given by the equation 
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Dt 


k£ 


3v„ 3v, 


- T, Tv~ - T , 


Dt km 8x £m 3x 

m m 


•k — ( u u. u ) 
3x v m k £' 


m 


3v 


3 v B i . 3v k 


Of A/ iN. x f iV 

T km dx + T £m 3x 2K ^ 3x T mn T n£ 
mm m 


3v, 


T T J +H. T . +H 0 T . + C e6. 


3x mn nk km m£ £m mk o k£ 
m 


(52) 


where the reader is referred to Haworth and Pope'*’'*' for the precise form of 
H. It is a simple matter to show that in a rotating framework (with angular 
velocity ft = ftk) , Eq. (52) is given by 


Dt M *1 . 3 . . 

Dt T km 3x “ T £m 3x ~ 3x ^ U m U k U £r 
m mm 


. . fv _ h T T 

+ T km Sx + T £m 3x ' 2K ^ 3x m nl 
mm m 


3v 


£ 


+ -sr — T t 1 )+H.t„+H„t 1 + C e<5 
dx mn nk km m£ £m mk o k£ 
m 


- B,ft[e,. n T. + £_, T n --^(£-1 T T „+£_„ T T,)] (53) 

1 1 3£m km 3km £m 2K ^ 3km mn n£ 3£m mn nk 


Haworth and Pope ^ proved that the inertial term in (53) vanishes for a 
two-dimensional turbulence and, hence, this Langevin model satisfies the 
principle of material frame-indifference in the two-dimensional limit unlike 
all other existing second-order closures. However, these authors also 
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showed that this model does not give rise to a Taylor-Proudman reorganization 
in a rapidly rotating framework. In fact, it can be shown that for a turbulent 
channel flow in a rapidly rotating framework, (53) yields a stress tensor 
of the form (45). While this is a substantial improvement over what is 
predicted by all other existing second-order closures, it is still not 
completely correct in that should vanish in the limit as Q -> °°. 

It will now be demonstrated that by making a small modification in the 
model for the tensor G proposed by Haworth and Pope 11 , a second-order closure 
can be obtained that gives rise to a Taylor-Proudman reorganization in rapidly 
rotating turbulent flows, thus, alleviating the problems discussed above. 

This proposed modification to (51) is as follows: 


G , 0 = -j - pj | W W t - K W W | 1/2 6, . 
k a YC 2 mn n P P m mn nm k* 

3v, . 3v, 

* t 2 * 6 ! 5 -5?: - TT B 1 — V * 

x, m 

where 

— 1 ^ V k ^ V £ 

W W - * 3^ " ^ 


(54) 


(55) 


and S 2 is a dimensionless constant (for B 2 = 0, this model reduces to that 
of Haworth and Pope). Since the new term containing the coefficient S 2 
vanishes in a two-dimensional turbulence, it follows that this model also 
satisfies the principle of material frame- indifference in the limit of 
two-dimensional turbulence. In a rotating framework with an angular 
velocity ^ = fik, this modified expression for G gives rise to a Langevin 


model of the form 
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° T )ct \ _ 3 . , 

Dt ~ T km 8x " T {jn 3x 9x ^ U m U k U £/ 

m m m 


3v £ 3v k 3 i 3v k 

3 l^km 3x + T £m 3x ^ 2K ^ 3x T mn T n£, 

mm m 


9v ji 

+ 3x T mn T nk^ + ^km T m£ + ^£m T mk + ^o^k#, 
m 


3 l^ £ 3£m T km + £ 3km T m& Tkf- £ 3km T mn T n& 


2B ? - 

+ e _ 0 t t , ) ] + — t-t l(w - e„ ft)( w 
3£m mn nk^ J ^1/ 2 lv mn 3mn np 

- e, ft)x - KTW - e, f2)(W - e_ 0,)\ l ^ 2 x, n 

3np J pm ^ mn 3mn v nm 3nm 1 k£ 


(56) 


If we consider the case of a rapidly rotating framework (i.e., the limit as 
n °°) , equation ( 56) yields the constraint 


3 1 [ £ 3£m T km 


e 3km T £m 


2K ^ £ 


3km T mn T n£ + £ 35,m T mn T nk 


)] 


26 2 

~p2 l £ 3mn £ 3np T pm 


£ 3mn £ 3nm ' 


1/2 


V 0 


( 57 ) 


A straightforward, although somewhat tedious calculation shows that (57) 
is satisfied if and only if the Reynolds stress tensor takes the two- 


dimensional form 
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XX 

V 

xy 

T yy 

0 

0 


0 - 
0 

0 - 


( 58 ) 


The key result that x must be zero can easily be seen by examining the 

4p Id 

t 

zz-component of (57) which reduces to 


2 h. 

K 1/2 



(59) 


It is thus clear that with just one minor modification, the Langevin 
second-order closure of Haworth and Pope yields a turbulence model that is 
consistent with the Taylor-Proudman theorem and is, thus, suitable for the 
description of rapidly rotating turbulent flows unlike all other existing 
turbulence models. 


No constraint is placed on x , x , and x 

xx xy yy 

for a two-dimensional turbulence. 


since Eq. (57) vanishes 



20 


4. Conclusion 

It has been proven that all existing second-order closure models 
are fundamentally inconsistent with the Navier-Stokes equations in a 
rapidly rotating framework. In particular, it was demonstrated that, 
for the test problem of fully-developed turbulent channel flow subjected 
to a rapid spanwise rotation, all existing second-order closures yield 
the spurious physical result of a vanishing Reynolds shear stress and are 
inconsistent with the Taylor-Proudman theorem. The type of results 
predicted for this problem by these second-order closures are, in fact, 
substantially more unphysical than those obtained from the simpler K-e 
or K-S, models of turbulence as indicated in Section 3. By adding one 
term to the second-order closure of Haworth and Pope obtained from a Langevin 
model, it was demonstrated that these problems can be overcome. This 
modified Langevin model appears to yield the first Reynolds stress closure 
that is generally consistent with the Navier-Stokes equations in a rapidly 
rotating framework. 

Future calculations should be conducted on rotating turbulent 
flows with this modified Langevin model. This is quite important since 
the closure model for G given in (54) is not the most general and, hence, 
there may be other terms which could have been overlooked. Unfortunately, 
such a computational study is quite involved and is beyond the scope of 
the present paper. Future research is also needed on the effect of 
rotations on the evolution of the joint probability density function since 
such evolution equations play an important role in all Langevin models. 

Few, if any, studies along these lines appear to have been reported. 
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With some additional refinements, turbulent closure models can be 
developed that apply to rotating flows for a wide range of rotation rates 
without the need for ad hoc empiricisms. Such a development could have 
a profound effect on the analysis of a wide variety of turbulent flows that 
are of interest to engineers and geophysicists. 
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FIGURE 1. 


Fully-Developed Turbulent Channel Flow Subjected to a 
Spanwise Rotation after Johnston, Halleen, and Lezius (1972). 
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